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We study the time evolution of correlation functions in 
closed quantum systems for nonequilibrium ensembles of ini- 
tial conditions. For a scalar quantum field theory we show 
that generic time-reversal invariant evolutions approach equi- 
librium at large times. The calculation provides a first princi- 
ples justification of Boltzmann's conjecture that the large- 
time behavior of isolated macroscopic systems can be de- 
scribed by thermal ensemble averages. 



I. INTRODUCTION 

One of the fundamental questions in statistical physics 
concerns the observation of macroscopic irreversible or 
dissipative processes arising from microscopic reversible 
dynamics. A classical approach to this question employs 
the separation of a system into observed and unobserved 
or environmental degrees of freedom and an averaging 
procedure over the latter Q). The reduced system may 
lose the detailed information about initial nonequilib- 
rium conditions — a necessary condition to reach ther- 
mal equilibrium. In contrast, in a closed system with a 
unitary time evolution no information is lost and it can- 
not exhibit dissipation or reach thermal equilibrium at a 
fundamental level. 

In this Letter we study the time evolution of closed 
quantum systems: for a real, scalar quantum field 
theory we show that generic time-reversal invariant 
evolutions of nonequilibrium quantum fields approach 
thermal equilibrium. We compute time-reversal in- 
variant evolution equations for correlation functions 
Tr {p(to) $(ii) . . . $(<„)}, i.e. expectation values of prod- 
ucts of Heisenbcrg fields $ for given initial density matrix 
p(to). We explicitly calculate the large-time limit numer- 
ically for a set of different initial conditions with given 
average energy density. We observe that the correlation 
functions approach the thermal distribution asymptoti- 
cally 



lim T>{p(to)$(tO 

|tO-*i|-»~ 



Mtn)}~Tr{e-' 3H <f>(t 1 )...<S>(tn)} 



without reaching equilibrium at accessible finite times. 
Typical evolutions of correlation functions for spatially 
homogeneous fields in 1+1 dimensions are shown in fig. 2. 
Our calculation provides for a quantum field theory a first 
principles justification of Boltzmann's conjecture: the 
large-time behavior of isolated macroscopic systems can 
effectively be described by equilibrium ensemble aver- 
ages. Apart from its fundamental meaning the treatment 
of closed quantum systems is of considerable practical 



importance. In many situations a clear separation into 
reduced system and environment is not obvious. This 
applies in particular to strongly coupled theories or sys- 
tems far away from equilibrium where a strict separa- 
tion of scales often does not exist. The nonequilibrium 
description of many-body quantum systems or quantum 
field theories has most diverse applications ranging from 
mesoscopic quantum devices to relativistic heavy-ion col- 
lision experiments. 

Nonequilibrium quantum field theory has attracted 
much interest in recent years. Nonpcrturbative ap- 
proaches based on mean field or large- N approximations 
as in [|]-Q or the use of evolution equations for gener- 
ating functionals || can give first principle insights into 
the dynamics of quantum fields. However, the large-time 
behavior of nonequilibrium quantum field theory posed 
an unresolved problem HU]. We study the time evolu- 
tion of correlation functions in a mean field approxima- 
tion, and in the Born collision approximation first sug- 
gested by Kadanoff and Baym for nonrelativistic real- 
time thermal Green's functions ||. Both descriptions 
lead to causal and manifestly time-reversal invariant evo- 
lution equations. While the mean field approximation 
fails to describe the long-time behavior j^-j^Q, we show 
that thermalization can be described in the collision ap- 
proximation without further assumptions. 



II. METHOD: LOOP EXPANSION OF THE 2PI 
EFFECTIVE ACTION 

We investigate a quantum field theory for a real, scalar 
field ip with classical action 

SM=Jd d+1 x(^d^{x)dMx)-^\x)-^\x)j (1) 

The time evolution of n-point correlation functions for a 
given initial density matrix p(to) can be obtained from 
the generating functional 

Z p [J, K] = Tr |p(t„) T c cxp % - (J J(x)$(x) 

K(x,y)$(x)<f>(y))\ (2) 

C,xy / ) 

by taking derivatives of Z with respect to the external 
sources J or K. In our notation J c = J d d x J c dxo 
and T c denotes time ordering along a closed time contour 
C starting at to and chosen to support </, K at times 
of interest lal. The introduction of the external bilocal 
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source term in (||) is used to construct the generating 
functional for 2PI Green's functions ||. The evaluation 
of the 2PI generating functional for a closed time path 
has been presented by Chou et al. 1 10 and by Calzetta & 
Hu jnj, and we refer the reader to their work for details. 
The path integral representation for Z p [J, K] 11 1C[| 

Z p [J,K] = Z[J,K- 1 ai,a 2 ,a 3 ,a i , . . .] = Jvtpexp ^(^M 
(J(x)+a 1 (x))ip(x) + - / (K(x ) y)+a 2 (x,y))ip(x)<p(y) 

C,x L JC,xy 
If 1 

+ ^t/ ai 3 (x,y,z)<p(x)<p(y)<p(z) + — 



3! JC,xyz 

<p(x)<p(y)ip(z)<p(w) 



4! 



a 4 {x,y,z,w) 



C,xyzw 



(3) 



exploits the fact that the information contained in the 
density matrix p(to) can alternatively be described by 
specifying all initial correlations. Here the source terms 
cti are non-vanishing at initial times t = to only. (An 
explicit representation of the sources in terms of the den- 
sity matrix is given in 0, section 6.) The 2PI ef- 
fective action is defined as the Legendre transform of 
Z[J, K\ ai, q;2, q;3, a4, , . . .] with respect to J and K M 



T[<j), G; «3, a4, ■ ■ ■] = —ih\nZ[J, K; 03, 04, . 
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(4) 

x)J{x) - 5 / U{x)4>{y) + HG(x, y)) K(x, y) . 

C,x Z JC,xy K ' 

Here we have absorbed a,\ into J and ai into K. The ef- 
fective action (Q) is parametrized by the macroscopic field 
<j)(x) = —ihS(\nZ)/5J(x) — (<p(x)) and the exact con- 
nected propagator KG(x,y) = (ip(x)tp(y)) — (ip(x))((p(y)) 
in the presence of the initial-time source terms 0:3, 04 etc. 
Solutions for <fi and G require 

6T[(f), G; a 3 , a 4 , . . .] 



5(f>(x) 
6r[<j), G; a 3 ,a 4 , 



-J{x) 



K(x,y)(f>(y) , 



C,x 



1 



hK(x,y) . 



(5) 



SG(x,y) 

where the RHS is zero for vanishing external sources J, K 
(t > to). Eq. (||) yields time evolution equations for the 
macroscopic fields (j> and G. We note that these equations 
are obtained from an action functional by a variational 
principle which guarantees a unitary time evolution. We 
restrict the discussion to a quartic initial density matrix, 
i.e. a.i = for i > 5, which is no approximation but 
constrains the initial state. We approximate the effec- 
tive action (Q) by a series expansion in orders of h or, 
equivalently, an expansion in the number of loops for 
two-particle irreducible graphs |||ll],[l(J . For simplicity 
we consider in this section the three-loop result for the 
effective action for a vanishing field expectation value <f>. 
Deviations from = can be included in a straightfor- 
ward way and results are described below. For the Z 2 - 
symmetric theory this implies a% = and the effective 
action T[G; a 4] = T[0, G; 0, 04] reads to order h 3 



r[G;a 4 }=- /[lnG- 1 ]^)-- 6 x , y (D x + m 2 )G(x,y) 



C.x 



9 / -*,v\ 

■'C,xy 



- -Q- / (^Sx, y 8x,z$x,w ~ a 4 (x,y,z,w))G(x,y)G(z,w) 

" JC,xyzw 

h 3 f 

~ 48 / (^x,yS x ,zS x ,w ~ ot4(x,y,z,w))L(x,y,z,w) (6) 

" J C,xyzw 

with the four-point function 



L(x,y,z,w) = -i \ (\5 x i t y'8 x i tZ >5 x > jW > 

JC.x'y' z'w' 

-a A {x , y', z\ w'))G(x, x')G{y, y')G(z, z')G(w, w') . (7) 



III. NONEQUILIBRIUM TIME-REVERSAL 
INVARIANT DYNAMICS 

The time evolution equation for the time ordered prop- 
agator 

G(x, y) = G> (x, y)Q(x - y ) + G> (y, x)Q{y - x ) (8) 

is obtained from (^) and the equal-time commutation re- 

lation i -£- a ( G >( x -.y) ~ G >(y^ x ))\x = yo = § x,y for vanish- 
ing external source K 

Xh 



U X G > (x, y) = -\ m + — G> (x, x) )G > (x, y) 



A/r 
IT 



L{x,y,y,y) 



(9) 



with the four-point function 



L(x,y,y,y) = -i\ / d z I I dzoG > (y,z)G > (x,z) 



+ dz G > (z,y)G > (x,z) - dz G> (z, y)G > (z, x) 
Jyo Jo 

+ L (x,y,y,y). (10) 



Here Lq is the initial-time four-point function which is 
determined by (Q) for given initial a 4 and G. Note 
that Gt,(x,y) = G > (y,x). The evolution for G > (x,y) 
is time-reflection invariant and causal since the "mem- 
ory integral" in ( |l0| ) only depends on times smaller than 
max(i ,i/o). 

The evolution equation for G> is a nonlinear, integral- 
differential equation which can be solved numerically. We 
use a standard lattice discretization for a spatial volume 
V with periodic boundary conditions and study the large 
volume limit to remove finite size effects. The time dis- 
cretization respects time-reversal invariance. Numerical 
results are calculated for 1+1 dimensions. All quantities 
will be expressed in units of appropriate powers of m with 
h = 1. We consider spatially homogeneous fields such 
that G> (x, y) = G> (xo, yo; x — y) and Fourier transform 
with respect to the spatial variables. The initial condi- 
tions are specified in terms of the momentum modes of 
the propagator G>(0, 0;p), its first derivative with re- 
spect to time G, and the four-point function Lq at initial 
time. 
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A. Two-loop (mean field) approximation 

The solution of (|9|) for the propagator to order 
h° corresponds to the free field solution and the 
propagator modes G>(i, 0;p) oscillate with frequency 
a/p 2 + m 2 /27r. The two-loop contribution to the effec- 
tive action (jfj) adds a time dependent mass shift Am(t) — 
IT I (2w) d C> Q) to the free field evolution equation. 
As a consequence the equation (||) for G becomes non- 
linear. The mass shift term in the evolution equation 
is the same for all Fourier modes G>(t, t';p) and each 
mode propagates collisionlessly with a time dependent 
effective mass. Fig. 1 shows a typical time evolution of 
the mass shift Am(t) to order h. As initial conditions we 
use the mean field thermal solution for G(0, 0; p) with 
inverse temperature (3 — 0.5 and G(0, 0; p) chosen to de- 
viate from the thermal solution. For Am(t) one observes 
an initial effective damping of oscillations in the mean 
field approximation termed dephasing The damp- 
ing is more efficient in higher dimensions Since 
the evolution equation is time-reversal invariant dissipa- 
tion is absent and the effective damping is due to the su- 
perposition of oscillatory functions with continuous fre- 
quency spectrum. This absence of dissipation can be 
easily demonstrated for finite systems with a discrete fre- 
quency spectrum. In this case the time-reversal invariant 
evolution equation, which is local in time to order h, leads 
to characteristic recurrence times after which the effec- 
tive damping is lost. Since we use a lattice discretization 
for a spatial volume V we have a discrete frequency spec- 
trum for finite volumes. In 1+1 dimensions we explicitly 
verify that the observed recurrence times for Am(t) scale 
with V or the number of degrees of freedom to infinity. 

The lower curve in fig. 1 shows a typical evolution of 
the macroscopic field 4>(t) in the mean field approxima- 
tion with nonzero initial 4>. The evolution of the field, but 
also the evolution of Fourier modes of correlation func- 
tions, are typically undamped in 1+1 dimensions. The 
mean field approximation does not describe the approach 
to a thermal equilibrium distribution [^| [^J . 
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FIG. 1. Mass shift (upper curve) and macroscopic field 
(lower curve) in mean field approximation. 



B. Three-loop (collision) approximation: 
thermalization 

A crucial improvement for the description of the long- 
time behavior comes from the three-loop contribution to 
the effective action ([|) . The resulting effective four-point 
function ( |Io| ) in the evolution equation for G> introduces 
scattering. In fig. 2 we show the time dependence of the 
equal-time propagator G> (i, t; p) for three Fourier modes 
p = 0, 3, 5 and three different initial conditions with the 
same energy density. For the solid line the initial con- 
ditions are close to a mean field thermal solution with 
(3 = 0.1, the initial mode distribution for the dashed and 
the dashed-dotted lines deviate more and more from ther- 
mal equilibrium. It is striking to observe that propagator 
modes with very different initial values but with the same 
momentum p approach the same large-time value. The 
asymptotic behavior of the two-point function modes are 
uniquely determined by the initial energy density which 
is characteristic for thermal equilibrium. 
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FIG. 2. Momentum modes of the two-point function for 
different initial conditions with same energy. 

We explicitly verify that energy is conserved and that 
the higher n-point functions approach asymptotic values 
uniquely determined by the energy. Note that all higher 
correlation functions are expressed in terms of G and 
their effective thermalization results from the behavior 
of the two-point function. We find that small deviations 
from = approach zero in the long-time limit for the 
initial conditions considered in fig. 2. 

The evolution of correlation function modes exhibits 
an effective damping of oscillations. The decrease of the 
maximum amplitude quickly approaches an exponential 
behavior, but the oscillations never damp out completely. 
This can be seen from the logarithmic plot of the prop- 
agator zero mode G>(i,£';0) in fig. 3. The upper curve 
shows |G>(i,0;0)| and exhibits a strong effective reduc- 
tion of correlations of the field at time t with the initial 
field. The lower curve presents the equal-time propa- 
gator G>(i,i;0), subtracted by its time average over 20 
time steps Gt and limited below to make it suitable for 
a logarithmic plot. 
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FIG. 3. Logarithmic plot for the propagator zero mode. 

We emphasize that time-reversal invariance is mani- 
fest, i.e. for each evolution towards the thermal solution 
there exists the reversed evolution away from equilib- 
rium. However, if one follows these reversed evolutions 
generically the system approaches equilibrium again after 
large enough times. We note that recurrence times can 
be infinite in the collision approximation even for small 
volumes V with discrete Fourier momentum modes. This 
is possible because the evolution equation for the ensem- 
ble average (||),(|l(J) is nonlocal in time. The memory ef- 
fects quickly lead to contributions which can be described 
by the superposition of oscillatory functions with a con- 
tinuous frequency spectrum. For the initial conditions 
employed in figs. 2 simulations for rather small volumes 
of order tens the correlation length are found to describe 
the large-time behavior well. 

The efficient damping observed in figs. 2 is due to a 
large coupling constant A = 10. The damping time r, 
defined as the time for which the envelope amplitude of 
G>(t, t; 0) is reduced by a factor e, is presented in fig. 4 
as a function of A. We find that the critical coupling for 
which the amplitudes are undamped is zero. The qualita- 
tive thermalization behavior does not differ in the weak 
and the strong coupling regime in the present approx- 
imation, which neglects higher-order scattering effects. 
These are relevant, in particular, for the computation of 
the bulk viscosity or the description of critical phenom- 

enaHS 

Of course, there are examples, like an energy eigenstate 
as initial condition, where thermalization can not occur. 
For few body systems like a set of coupled anharmonic 
oscillators it is known that a transition to non-ergodic 
behavior can occur for small but nonzero couplings [fl3[ . 
For a finite volume V and ultraviolet cutoff for the mo- 
mentum modes our model can be mapped onto a finite set 
of coupled anharmonic oscillators. We indeed find that 
for small enough volumes and couplings below a certain 
critical coupling A c the long-time evolution can be un- 
damped. The values for the critical coupling strongly 
depend on the initial conditions and the volume or num- 
ber of degrees of freedom. 
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FIG. 4. Damping time r for propagator modes with mo- 
menta |p| = 0, 3, 5 (from top) as a function of the coupling A. 

To verify the equivalence of the large-time limit of 
nonequilibrium correlation functions with the description 
in terms of a thermal distribution, we directly compute 
equilibrium correlation functions for p = exp(—(3H). 
The computation of the quantum-statistical partition 
function is done in the standard imaginary-time formal- 
ism, equivalent to evaluating (Q) for an imaginary time 
path C running from zero to —i/3. From the stationary 
condition for the corresponding imaginary-time three- 
loop effective action one obtains the gap equation for the 
propagator Fourier modes 

G-/ ( W , , p) = wf + p 2 + m 2 + J2j (^^K, q) 

G E {ui - uj n - Lu m ,p - q - k) . (11) 

Here oji — 2ttI/(3, 1 6 Z and the relation between the Eu- 
clidean propagator Ge and the equal-time propagator G 
is given by TJ2i Ge(^i,p) = G(t — t;p). We solve eq. 
( pd| ) by iterating the RHS, starting from the mean field 
solution, until convergence is obtained. We find that 
the large-time result for the equal-time correlator from 
the real-time dynamics and from the partition function 
agrees in the continuum limit. Their correspondence pro- 
vides a first principles justification of Boltzmann's con- 
jecture. 
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